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Abstract. We calculate the masses of the resonances D* (2317) and D s i(2460) as well as their 
bottom partners as bound states of a kaon and a and B^'-meson, respectively, in unitarized 

chiral perturbation theory at next-to-leading order. After fixing the parameters in the D^o(2317) 
channel, the calculated mass for the D s i(2460) is found in excellent agreement with experiment. 
The masses for the analogous states with a bottom quark are predicted to be Mb* = (5696 ± 
40) MeV and Mb b1 = (5742 ± 40) MeV in reasonable agreement with previous analyses. In 
particular, we predict Mb s1 —Mb* = 46 ± 1 MeV. We also explore the dependence of the states 
on the pion and kaon masses. We argue that the kaon mass dependence of a kaonic bound state 
should be almost linear with slope about unity. Such a dependence is specific to the assumed 
molecular nature of the states. We suggest to extract the kaon mass dependence of these states 
from lattice QCD calculations. 

PACS. 12.39.Fe Chiral Lagrangians - 13. 75. Lb Meson-meson interactions - 14. 40. Lb Charmed 
mesons 



1 Introduction 

The discovery of two narrow resonances with open 
charm in 2003 opened a new chapter in hadronic spec- 
troscopy. The D* (2317) was discovered by the BaBar 
collaboration [T] and shortly after the Z? s i(2460) by 
the CLEO collaboration g]El Both were confirmed 
later. 

These resonances were not in agreement with the 
predictions of constituent quark models, for instance 
the Godfrey- Isgur model [3], that on the other hand 
were very successful for all other open charm states 
discovered so far. The masses of the Z?* (2317) and of 
the D s i(2460) are approximately 160 MeV and 70 MeV, 
respectively, below the predictions in Ref. [3]. In ad- 
dition, 

■^n 3 i(2460)-^n* (2317) — M D *-M D . (1) 

This equality is difficult to explain in a conventional 
quark model. It can be explained in a parity dou- 
bling model [415] . This result is preserved at one-loop 

1 BaBar already observed a signal of the D s i (2460) in [I] . 



level [5]. However, it is not clear if the parity doubling 
is a consequence of quantum chromodynamics (QCD). 
The equality of mass splittings appears naturally, if 
both states are assumed to be molecular states of a 
kaon and aflora D* [7], since spin-dependent interac- 
tions for a heavy quark are suppressed by two powers 
of Aqcd/uiq (see discussion in Section [2]), with Aqqd 
the soft QCD scale and rriQ the mass of the heavy 
quark. Soon after the discovery, based on phenomeno- 
logical calculations, Barnes and Close [5] and van Bcv- 
eren and Rupp [9] pointed out the possible molecular 
nature of these resonances: they described them as 
bound states of D^*'K mesons. For reviews of possi- 
ble explanations of these states, see Refs. [10111] . The 
purely phenomenological calculations were shortly af- 
ter improved by Refs. |12ll3ll4ll5j . where variants of 
unitarized chiral perturbation theory were used, and 
in addition open bottom molecules were also studied. 
In Refs. [16117] both resonances emerged as molecules 
from Weinberg- Tomozawa type interactions using flavor- 
SU(4) constraints. 

Clearly it is important to find observables that al- 
low one to identify the nature of the mentioned states 
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model-independently either from data or from a com- 
parison to lattice QCD results. Up to now a few such 
quantities were pointed out, namely the value of the 
hadronic width of the resonances [18119120] — which 
in a molecular assignment turns out to be significantly 
larger than that for a quark state [21] — and the 
value of the KD scattering length in the resonance 
channel [32] — which gets large and negative, if the 
states are molecules [23124] . In this work we add an 
additional quantity, namely the quark mass and es- 
pecially the kaon mass dependence of the resonance 
masses. We demonstrate that this dependence is qua- 
litatively different for quark states and for hadronic 
molecules and thus, a study of this dependence using 
lattice QCD is of high importance. 

Our study is based on Refs. [20122) . but extended 
in that now D- and D*-mesons are both included on 
equal footing. This allows us to also investigate the 
implications of heavy quark spin symmetry. 

The paper is organized as follows. In Section [2] 
we give the Lagrangian for and display the calcula- 
tion of Goldstone boson scattering off D- and D*- 
mesons as well as the unitarization scheme used to 
generate the £>* (2317) and the D s i(2460) as hadronic 
molecules. In Section [3] the coupled channel problem 
in the strangeness 5 = 1 isoscalar channel is solved 
and the D* (2317) and the D sl (2A60) are generated 
as DK and D*K molecules. We discuss the extra- 
polation of their masses to unphysical values for the 
pion and kaon mass, respectively. A brief summary 
and outlook is given in Section [4] The appendix con- 
tains some details on the derivation of the pertinent 
scattering amplitudes. 



2 Scattering amplitudes 

The leading order (LO) Lagrangian in Heavy Meson 
Chiral Perturbation Theory is given by the kinetic en- 
ergy of the heavy fields, the coupling of the heavy 
fields to pions and the mass splitting of the heavy 
mesons [25126127] 



where 



+—Tr[H a a^H a a^] 
with the heavy fields defined as p 



(2) 



H 



1 



[Y + iPls], H = 7°Hh°, 



where Latin indices denote flavor and the trace is 
taken over the Gamma matrices. Here, the heavy fields 
are P = (D°,D+,D+) and = {Df, D*+, D* s +), 
and = (1, 0) + 0(p/M) is the heavy mesons velo- 
city. The covariant derivative is 



U = exp 



U 



(4) 



with F the pion decay constant in the chiral limit. Up 
to the order we are working, F can be replaced by the 
physical pion decay constant F n = 92.4 MeV. As we 
will see, it is basically F that provides the strength of 
the interaction. 

The Goldstone boson fields are collected in the ma- 
trix <t> with 



7T 

K 



r\ 7r + K + 



The axial- vector current is defined as 



,t 



(5) 



(6) 



The low-energy constant (LEC) g v can be fixed from 
the strong decay of the D* + . Taking the values from 
the Particle Data Group (PDG) [29], we find g n = 
0.30 ±0.03. 

The spin-symmetry breaking term gives rise to a 
mass difference 



A = Ah 



M P = -i 



(7) 



with rriQ the pertinent heavy quark mass. In this work 
we use physical values for the meson masses and thus 
the mentioned term is considered automatically. Note, 
since M Da * - M Ds jt M D * - M D (c.f. Eqs. dTUJl), in 
this way we also include an effect of simultaneous spin 
symmetry and SU(3) violation, which is formally of 
next-to-next-to leading order (N 2 LO). We come back 
to the quantitative role of this sublcading effect below. 
The Lagrangian at next-to-leading order (NLO) is 

■5-nlo = + r R[H a (-h (x+)bb + h 2 {u fl u ,1 ) bb 
+hi(-iv ll )(u^u v ) b h(iv u ))H a ] 
+Tr[H a {h 1 (x+) ba + h 3 (u^) ba )H b ] 
+h 6 Tr[S a (-iv li ){u' i , u»} ba (iv u )H b ], (8) 

where we have introduced explicit chiral symmetry 
breaking terms due to non- vanishing light quark masses 
using 



X± = " t xu t ± ux^u , 



X± = X+ 



+ )aa i 



(9) 

\qq\0)\/F 2 and 



i 



with x = 2B-diag(m„,md,m s ), B = 
(x+)aa is the trace of x+- 

Following Ref. [T^j] we neglect the terms ho, h 2 
and hi since they are suppressed in the large- Nc limit 
+ ud^v)) (3) because of the second flavor trace. The LEC hi can be 
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Fig. 1. Diagrams contributing at LO and NLO to Goldstone boson scattering off D- and D*-mesons and the pertinent 
kinematics. Dashed lines denote Goldstone bosons, solid lines charmed pseudoscalar mesons and solid double lines 
charmed vector mesons, in order. 



determined using the relation B(m s — m) — \^j- k o 
M^ + )/2— M^ as done in [13]. We use the value hi = 
0.42 ±0.00. 

Since we are interested in the masses of the reso- 
nances, isospin breaking can be neglected. Hence, we 
take averaged values for charged and neutral particle 
masses, and the pion decay constant from the PDG, 
as listed in the following [29] : 



which also preserves spin symmetry. Its S-wave pro- 
jection is 



vr =v 



yjMihh — 



-C35 — — h 35 E 2 E4 



(13) 



M D = 1867 MeV , 
M D * = 2008 MeV , 
M w = 138 MeV , 
M v = 548 MeV . 



M Ds = 1968 MeV , 
M D . = 2112 MeV , 
M K = 496 MeV , 



(10) 



We will calculate the scattering of the Goldstone 
bosons off the pseudoscalar D-mesons (0 + channel) as 
well as off vector D* -mesons (1 + channel). The dia- 
grams at the order we are working are shown in Fig. [TJ 
The diagrams are evaluated in the isospin basis, the 
relation between the isospin basis and the particle ba- 
sis can be found in the Appendix of Ref. [35] . 

Evaluating the diagrams for the leading order con- 
tact interaction we find that the contributions for the 
+ channel and the 1 + channel are the same up to 
the different masses for the pseudoscalar and vector 
charmed mesons, as expected from heavy quark spin 
symmetry: 



1 



C y/M 1 M 3 ^{E 2 



£4). (11) 



The constants Co for the various channels are listed 
in Tab. [T] Further, Mi and M 3 are the masses of the 
in-coming and out-going charmed meson, respectively, 
cf. Fig. [T] and E 2 and E4 are the energies of the in- 
coming and out-going light mesons. 

For the NLO contact interactions we again find 
the same contributions except for the different pseu- 
doscalar and vector charmed meson masses: 



V 0+ 



NLO 



v 1+ 



+C 35 {h 3 p 2 -p4 + 2h 5 E 2 E 4 ) 



(12) 



with /135 = h 3 + 2/15 = I135/ME). This is the only free 
parameter at this order. This is fitted to the mass of 
the D* (2317). We get h 35 = 0.35 for the dimension- 
less parameter, which is of natural size. 

The s-channcl exchange needs a P-wave interac- 
tion and thus does not need to be considered. The only 
contribution from exchange diagrams that does not 
have an vanishing S-wave projection is the w-channcl 
exchange of a charmed pseudoscalar meson in the 1 + 
channel: 



,2 r 



Ex 



— Cu 



3F 2 



Mi 



|P| 2 + 



E 2 

TT 3 



,'|2 



(y ■ k) 



VM1M3 



(14) 



where the constant C u is also given in Tab. [T] and p 
and p' are the three-momenta of the in-coming and 
out-going mesons in the center-of-mass frame and the 
energy transfer (v ■ k) is defined in Fig. [1] Observe, 
since k, p and p' are counted of order of the Ei, this 
potential appears formally at NLO. Since it acts in the 
1 + channel, while there is no counter part in the + 
channel, it thus provides formally the leading source of 
spin symmetry violation. However, in practice its con- 
tribution turns out to be very small, due to three rea- 
sons: first of all g% = 0.1, second, |p| ~ y/Mxe, with 



e ee M d ,+M k -M Db 



(2460) 



M D + M K -M 



if- JW ZT (2317) 



being the binding energy, which is significantly smaller 
than Mx and, most importantly, it does not operate 
in the KD* -> KD* channel (in that channel C u = 
— cf. Tab. [T]). In total it gives a negligible contribu- 
tion. 

For SU(3) calculations, the uncertainty from the 
chiral expansion up to NLO is & (M K /A^. ) , where A x w 
^ttF^ sets the hard scale in the chiral expansion, com- 
pared with the LO contributions. So we use 



V 



(15) 
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where V s is the S-wave projection of the full NLO am- 
plitude calculated here for the + and 1 + channel, 
respectively. A more detailed discussion of the scat- 
tering amplitudes is given in the Appendix. 

Dynamical generation of bound states from a the- 
ory without bound states as fundamental fields is a 
non-perturbative phenomenon. It cannot be provided 
by any finite perturbative expansion in the momen- 
tum. Thus we have to unitarize, i.e. resum, the ampli- 
tude we had obtained so far. Following Ref. [3D], the 
unitarized amplitude can be written as 



T( s ) = V ± (s){l-G(s)-V ± ( S )}- 



(16) 



where G(s) is the diagonal matrix with non- vanishing 
elements being loop integrals of the relevant chan- 
nels — see Eq. (fT7|> below. The bound state mass is 
found as a pole of the analytically continued unita- 
rized scattering matrix T . To be specific, the pole of a 
bound state, as is the case here, is located on the real 
axis below threshold on the first Ricmann sheet of the 
complex energy plane (note that because we neglect 
isospin-brcaking, these bound states do not acquire a 
width). In our approach, both the D* (2317) and the 
-D s i(2460) appear as poles in the S = 1 isoscalar 1 = 
channel. Since they couple predominantly to DK and 
D*K, we interpret the £>* (2317) and D s i(2460) as 
DK and D* K bound states, respectively. 

At this point we have to make some comments on 
the loop matrix. In previous works on this subject 
[16117122] . the relativistic two-particle loop function 



d 4 q 



1 



1 



(2tt) 4 q 2 -M 2 +ie {q-P) 2 -M 2 



(17) 



was used in dimensional regularization with a subtrac- 
tion constant a(/i) which absorbs the scale dependence 
of the integral, where M$ and Mn are the masses of 
the light meson and charmed meson involved in the 
loop and P denotes the total momentum. However, 
this relativistic loop function violates spin symmetry 
very strongly. Particularly, using the parameters fitted 
to the mass of the Z?* (2317), we find the mass of the 
£> s i(2460) at 2477 MeV. From a field theoretical point 
of view, using a relativistic propagator for a heavy 
meson and using the Lagrangian from heavy quark 
expansion simultaneously is inconsistent. There is a 
well-known problem in relativistic baryon chiral per- 
turbation theory for pion-nucleon scattering, which is 
closely related to the case here. The relativistic nu- 
cleon propagator in a loop explicitly violates power 
counting, see [51] . In our case this effect appears only 
logarithmically as log (M|,/^ 2 ). However, it induces a 
violation of spin symmetry in the interaction already 
at NLO and thus violates the power counting. 

There arc different ways to deal with this. In [T^] a 
subtraction scheme is used that identifies the scale /j, 
with the mass of the D- and D*-mesons, respectively, 
and thus the spin symmetry violating terms disappear. 



In Refs. [16117] simply different subtraction constants 
are used for the two channels. In this work we use a 
static propagator for the heavy meson, as this is con- 
sistent with the Lagrangian in the heavy quark expan- 
sion. This also allows us to use the same subtraction 
constant and scale of regularization for both channels. 
The heavy meson propagator for the vector D*-meson 
then takes the form 



{P»",l} 



{P^,!} 



Ml 



it 2M D (vk-A 



te 



(18) 



where P Miy is a projector for spin-1. Dropping A in 
the denominator and using the one in the numerator 
of Eq. ([T5]) amounts to the propagator for the pseu- 
doscalar £)-meson. 

With this the integral / defined in Eq. ([IT]) turns 
into 



I 



1 



16n 2 M D 



(P°-M D ) 



( M T 

a(ju)+log I — f 



-2 y /(P a -M D ) 2 -M 2 cosh- 



_! (P°-M D 



Ma, 



-27Tl|p c 



o 



M 



D 



(19) 



with a(/i) a subtraction constant, \i the scale of dimen- 
sional regularization, which is fixed to the averaged 

mass of D and D* , and |p cms | = yJ(Po - M D ) 2 - Mj. 

Note, the use of the modified loop function of Eq. (fl^j) 
leads only to tiny changes from relativistic effects in 
the observables discussed in Refs. [20122] . The sub- 
traction constant a(n) is fitted to the mass of the 
D* (2317). We find a{fi = 1936 MeV) = -3.034. Us- 
ing the same parameters for the axial-vector channel 
we find the pole at 



M 



D s i(2460) 



(2459.6 ±0.6 ±1.8) MeV 



(20) 



where the first uncertainty stems from the experimen- 
tal uncertainty for the mass of the D* (2317), which 
was used for fitting the parameters, and the second 
4% uncertainty estimates the higher orders from the 
heavy quark expansion — since there is no significant 
spin symmetry breaking term in the scattering ampli- 
tudes at NLO, spin symmetry breaking effect should 
appear at ^([^qcd/toc] 2 )- Hence the second uncer- 
tainty can be estimated by (v1qcd/to c ) 2 e. In this con- 
text it is interesting to note that the spin-symmetry 
violation induced by Mb s * — Md s ^ Md^ — Mjj, which 
is of N 2 LO, indeed contributes to Md s1 ( 2 46o) less than 
1 MeV, consistent with the uncertainty estimate. 

The result in Eq. ([2"0]) is in perfect agreement with 
the experimental result 



A/ cxp 

JW D 3 i(2460) 



(2459.5 ±0.6) MeV . 



To the order we are working instead of F w we could 
as well have used Fk = 1.18F ff . However, this change 
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Table 1. Coefficients for the amplitudes for all possible channels with total strangeness 5* = f and total isospin 1 = 0. 
While Co, Ci, and C35 act in both channels + and 1 + , while C u acts solely in f + . 



Channel 


Co 


Ci 


C35 




DK -> DK 
D a -q -> D a r) 
D a rj -» DK 


-2 


-Vs 


-4M| 

-2(2M 2 - Ml) 


2 

4 
3 

V3 
3 




2 
3 

_ ^/T 



Table 2. Comparison of our predictions of the masses 
of the BK and B* K bound states with those in 
Refs. [12I14I15J . All masses are given in MeV. The un- 
certainties given in the first column originate from the 
residual scale dependence and an estimate of higher or- 
der effects, respectively. 





This paper R e f. |12| Refs. 1141151 


Mb- 
Mb b1 


5696 ± 20 ± 30 5643 5725 ± 39 
5742 ± 20 ± 30 5690 5778 ± 7 



docs not have any impact on the mass calculated for 
the D s i(2460) once the subtraction constant is re- 
adjusted to reproduce the mass for the scalar state. 

The formalism used here can easily be applied to 
the meson sector with open bottom as well. In order 
to do so, we just have to replace the D-mesons with 
the corresponding 7?-mesons. Their masses are |29j : 

M s = 5279 MeV, M Ba = 5366 MeV, 

M s * = 5325 MeV, M s . = 5415 MeV. (21) 

Heavy quark flavor symmetry implies that all the other 
parameters stay the same up to higher order correc- 
tions._This allows us to predict the masses of the BK 
and B*K bound states in the (S = 1,7 = 0) channel. 
The predicted masses are given in Table [5J together 
with a comparison with previous calculations |12ll4ll5j 
Our results are in reasonable agreement with the pre- 
vious ones. In our calculation we have two sources of 
uncertainties, both shown explicitly in the table. One 
stems from higher order effects, which may be ob- 
tained from multiplying the binding energy by 20%, 
estimated as ©(/Iqcd/toc) since heavy flavor symme- 
try was used to relate the LECs in charm and bot- 
tom sectors. The other one originates from the in- 
trinsic scale dependence of the result: To come to 
a fully rcnormalization group invariant amplitude, a 
complete one-loop calculation would have been nec- 
essary for the transition amplitude. This, however, is 
beyond the scope of this paper. Thus, when connecting 
the charm to the bottom sector, a residual scale de- 
pendence remains — to quantify it we varied the scale 
parameter fi (see Eq. (|19[1 ) from the averaged mass of 

2 The B^ mesons in the calculations of Ref. [T5] should 
be understood as B'*' mesons which contain a b quark 
rather than b. 



D and D* to that of B and B* , while keeping a(/x) 
fixed. Combining the two uncertainties in quadrature 
gives a total uncertainty 40 MeV which is about 1 
% for the masses. When again switching from TV to 
Fk as the strength parameter of the interaction (c.f. 
Eq. Q and sentences below), the predicted masses 
change by 8 MeV only, well consistent with our uncer- 
tainty estimates. 

It should be stressed that the uncertainties quoted 
for Mb* and M~b s1 are highly correlated. As already 
explained for the charmed system, within the molec- 
ular scenario the relation 

M B Mir M B , -M B (22) 

should hold up to corrections of 6 \[Aqc d / 'mb} 2 ) ~ 
c.f. discussion below Eqs. (fl"4"|) and (f2"U|). Thus we pre- 
dict 

M Bsl -M B ; = 46 ± 0.4 ± 1 MeV , (23) 

where the first uncertainty comes from the current ex- 
perimental uncertainty in Mb-— Mb and the second 
from the estimated spin breaking effects in the for- 
mation of the molecule. Clearly, all mass differences 
deduced from Tab. [2] are consistent with this value. 

M Bsl and M B * o have not been measured so far. 
Note that their existence in the deduced mass range 
and, especially, with the mass splitting of Eq. (|23|). 
is a crucial and highly non-trivial test for the theory 
pbresented and especially for the molecular nature of 
ooth states. 



3 Chiral extrapolations 

To test the nature of the resonances besides further 
experimental data we can also compare our results to 
lattice calculations (for a corresponding study in the 
light meson sector see Ref. [35] ). To do so, we extend 
the calculations from the physical world to unphysi- 
cal quark masses which are frequently used in lattice 
calculations. Varying the light quark masses is equiv- 
alent to varying the pion mass. Although the physical 
strange quark mass is nowadays routinely used in lat- 
tice calculations, we emphasize that by varying the 
strange quark mass, or equivalently varying the kaon 
mass, one can learn a lot about the nature ol some 
hadrons, as will be discussed below. So we have to 
express all results in terms ol pion and kaon masses, 
respectively In the following chapter only the charmed 
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mesons are discussed explicitly, however, analogous ar- 
guments apply to their bottom partners. 

In order to proceed we need to assume that the 
subtraction constant a(/i) does not depend on the light 
quark masses. We stress, however, that even allowing 
for a quark mass dependence of a(fi) would not change 
the general features of the results, but might slightly 
enhance the uncertainties. 



3.1 Pion mass dependence 

Lattice QCD calculations are often performed at larger 
quark masses than realized in nature — for the KD 
system of interest here exploratory lattice studies are 
presented in Rcf. [33] , In addition, as we will argue 
here, the quark mass dependence of a state contains 
important information on its nature. Varying u and d 
quark masses can be expressed by varying the pion 
mass. In [34] . the charmed meson masses were ex- 
panded up to one loop order in the chiral expansion. 
Nevertheless for our purposes the expansion up to 
0(Ml) is sufficient. 

Using the Lagrangian given in Eq. ((5J) we find the 
NLO correction to the charmed meson masses to be 

5M 2 D{ , } = AhxBm , SM^ M = AhxBm s , (24) 

with rh = (m u + m c i)/2 the average light quark mass. 
When studying the pion mass dependence, we consider 
the physical value for the strange quark mass here, 
thus we can use the physical mass for the D s . Using 
Ml = 2Bm yields, up to 0(M%), 

M 2 

M D w=M DM +hi^-*-, (25) 

o 

where M D (*) is the charmed meson mass in the SU(2) 
chiral limit m u = rrid = with m s fixed at its phys- 
ical value. For the kaon and the eta mass we can 
find similar expressions by using M\ = B(rh + m s ). 

2 2 

M K = Bm s and M 2 = B|(m + 2m s ), M v = 5§m 5 
from the LO chiral expansion: 

M 2 o M 2 

M K = M K + -j 2 -. M v = M v + (26) 
AM k 6M V 

In Fig. [5] we show the mass of both resonances, 
£>* (2317) and D s i(2460), as well as their binding en- 
ergies as a function of the pion mass. Note that the ob- 
served rather strong pion mass dependence is specific 
for a molecular state. The corresponding dependence 
for a quark state should be a lot weaker. To under- 
stand this, one notices that a pure cs state does not 
contain any u, u, d 1 d quarks. The leading term con- 
taining these light quarks is l/N c suppressed. Thus, 
for the quark state, the pion mass dependence can 



only enter via D^K loops — as in case of the molec- 
ular state. However, their contribution should be a 
lot smaller for the quark state than for the molecular 
state. To see this, observe that the loop contributions 
from a particular meson pair is proportional to the 
squared coupling of that meson pair to the resonance. 
As shown in Refs. [23124) . this coupling is proportional 
to the probability to find the molecular state in the 
physical state. Thus, the pion mass dependence for 
£>* (2317) and £> s i(2460) should be maximal if both 
are pure molecules. This case is depicted in Fig. [5] On 
the other hand it should be very weak for the admit- 
tedly unrealistic scenario of a pure quark state. Note 
that the mentioned relation between the coupling and 
the structure of the state holds only in leading order in 
an expansion in \JMk^I ft (see Ref. [35] for a detailed 
discussion), with the inverse range of forces (3 ~ m p . 
For the scalar and axial vector open charm mesons 
this gives an uncertainty of the order of 20%. As a re- 
sult of the larger binding, for the bottom analogs we 
even find 30%. In addition in the analysis we neglected 
terms of 0(l/N c ). Thus, from this kind of analysis 
at most statements like 'the state is predominantly 
molecular/compact' arc possible. 

Furthermore we notice that the plots show an al- 
most identical behavior in the scalar and the axial- 
vector channel. So the spin symmetry breaking effects 
are only very weak here. We see that the binding en- 
ergy in both cases varies from about 40 MeV to about 
80 MeV. 



3.2 Kaon mass dependence 

Before going into details of calculations, let us make 
some general statements about the Mk~ dependence 
of the mass of a bound state of a kaon and some other 
hadron. The mass of such a kaonic bound state is given 

by 

M = M K + M h - e, (27) 

where Mh is the mass of the other hadron, and e de- 
notes the binding energy. Although both Mh and e 
have some kaon mass dependence, it is expected to 
be a lot weaker than that of the kaon itself. Thus, 
the important implication of this simple formula is 
that the leading kaon mass dependence of a kaon- 
hadron bound state is linear, and the slope is unity. 
The only exception to this argument is if the other 
hadron is also a kaon or anti-kaon. In this case, the 
leading kaon mass dependence is still linear but with 
the slope being changed to two. Hence, as for the D K 
and D*K bound states discussed here, we expect that 
their masses arc linear in the kaon mass, and the slope 
is approximately one. As we will see, our explicit cal- 
culations confirm this expectation. 

3 The /o(980) was proposed to be such a KK bound 
state |36I24] . 
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M„ [MeV] Mn [MeV] 




M n [MeV] Mn [MeV] 

Fig. 2. The masses and binding energies of the D* (2317) and the D s i(2460) as a function of the pion mass. 



For calculating the kaon mass dependence, we use 
physical masses for the pion and the charmed mesons 
without strangeness. To express the strange quark 
mass dependence in terms of kaon masses, we write 

M 2 K = B (m, + m) , M K = Brh = , (28) 



where M K is the mass of the kaon in the limit m s = 0. 
The charmed strange meson mass is 



M D w = M d h + 2hi 



M% Ml 
m— 



M D^ 



M D" 



(29) 



where M D (,) is the charmed strange meson mass in 



the limit m s = 0. Finally the eta mass is given to this 
order by the Gell-Mann-Okubo relation 



M* = B [ -m s 



m 



= -mI-\mI. 



(30) 



4 Certainly, lattice simulations use unphysical values for 
up and down quarks. However, the conclusion of this sub- 
section will not be affected by the value of pion mass. For 
definiteness, we choose to use physical masses for the pion 
and non-strange charmed mesons. 



For analyzing the kaon mass dependence we calcu- 
late the masses of the L>* (2317) and £> s i(2460) and 
their binding energies as before. The corresponding re- 
sults arc shown in Fig. [3] Furthermore we also show 
the D^K threshold in dependence on the kaon mass. 

For the Z?* (2317), we find that the kaon mass de- 
pendence of the mass is almost perfectly linear, espe- 
cially in the region M K = 300 - 500 MeV. For higher 
strange quark masses the chiral expansion is no longer 
valid and the results become less reliable. The kaon 
mass dependence of the mass and the DK threshold 
have a similar slope. The slope of the threshold is one, 
since the D-meson mass is independent of Mr- to the 
order we are working, and the slope of M D * (2317) is 
about 0.85. This finding is perfectly consistent with 
the expectation formulated in the beginning of this 
subsection. The deviation of the slope from one can 
be understood from the kaon mass dependence of the 
binding energy (see right column of Fig. [3]) and from 
the effects of the coupling to the D s ij channel. The 
situation for the -D s i(2460) is similar. 

On the contrary, if L»; o (2317) and L» s i(2460) arc 
assumed to be quark states, their masses necessarily 
depend quadratically on the kaon mass — analogous 
to Eq. (|2"9")l . Thus, an extraction of the kaon mass de- 
pendence of the masses of D* (2317) and D s i(2460) 
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Fig. 3. The blue band shows the masses and binding ener] 
mass, the red dashed line shows the threshold. 

from lattice data is of high interest to pin down the 
nature of these states. 



4 Summary 

In this paper we used unitarized heavy meson chi- 
ral perturbation theory to calculate the masses of the 
resonances D* (2317) and £> s i(2460) as well as their 
bottom analogs B* and B s i. The results for the latter 
turned out to be consistent with previous analyses. Es- 
pecially we predict Mb s1 —Mb* = 46 ±1 MeV , where 
we added the uncertainties of Eq. ([2"3"]) in quadrature. 
In addition, calculations were performed with varying 
pion and kaon masses. Within the molecular picture 
used here, we found a linear dependence of the res- 
onance masses on the kaon mass with the slope be- 
ing approximately unity and a strong quadratic de- 
pendence on the pion mass. We argued that, if both 
states are assumed to be elementary quark states, the 
mentioned dependences would be qualitatively differ- 
ent: while the pion mass dependence should still be 
quadratic, it should be a lot weaker, the kaon mass 
dependence should come out quadratic as well. Thus, 
especially an investigation of the kaon mass depen- 
dence of the £>* (2317) and D s i(2460) or their bottom 




300 350 400 450 500 550 600 



M K [MeV] 




300 350 400 450 500 550 600 



M K [MeV] 

;ies of D* (2?>Y7) and D s i(2460) in dependence of the kaon 



analogs within lattice QCD would be very interesting 
in order to determine the nature of these resonances. 
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A Amplitudes 

Here we give all the contributing amplitudes and their 
S-wave projections. In all the following expressions, 
Pi(3) and P2(4) denote the momentum of the in-coming 
(out-going) heavy meson and light meson, respectively, 
the heavy meson mass is given by Mu^\, and the en- 
ergy of the light meson is given by i?2(4)- The lead- 
ing order contact interactions are purely S-wave by 



Cleven et al.: Light meson mass dependence of the positive parity heavy-strange mesons 



9 



construction and identical up to the masses for both 
channels: 



1 



T L ° = ^Mjj 3 C — (E 2 + E A ) 



T i+ 



VMjf 3 (e ■ e*)C ^(E 2 + E 4 ) 



= ^M 1 M 3 C Q ^(E 2 +E 4 ), 

where e is the polarization vector of the D* 
The NLO contact interaction reads: 



(31) 



(32) 



rpO + 

J NLO 



y/MiM 3 



Ci — =hi 
3F 2 



+C 35 — (h 3 p 2 ■ p 4 + 2h 5 E 2 E 4 ) 



(33) 



J NLO 



+C 35 ^(h 3 p 2 ■ pa + 2h 5 E 2 E 4 ] 



Ci — -h x 

3F 2 



-C 35 j^(h 3 p 2 ■ p 4 + 2h 5 E 2 E 4 ) 



(34) 



The matrix element for the exchange of charmed 
vector mesons in the 1 + channel is zero at this order: 



1+ 4gl 1 

-Vector-Ex F 2 2vk 
= 0, 



(e • pe • p e • ep • p ) 



(35) 



where e is the spatial component of the polarization 
vector of the D* , and p^'' denotes the three-momentum 
of the in-coming (out-going) heavy meson in the center- 
of-mass frame and k is the residual momentum of the 
exchanged heavy meson (c.f. Fig. [1]). 

The exchange of charmed vector mesons in the + 
channel gives a non- vanishing contribution: 



- Vector — Ex 



-8C u y/hhNh 



ag p • p 

F 2 2v ■ k 



(36) 



But since it is purely P-wave, the S-wave projection 
vanishes. Finally we consider the exchange of charmed 
pseudoscalar mesons in the 1 + channel: 



T^s-Ex = -Cu^(e* •p 2 )(e-P4)^;v / ^TM3 



-c u 



M 

F 2 



M 4 M 3 
1 



^.| p |2 cos 2^ + g |prcos . 

IpIIp'I + IpIIp'I ) cos6»cos0' 



2(v ■ k) 



VM1M3, 



(37) 



where 9 (Q 1 ) is the angle between p (p') and the quan- 
tization axis used for the polarization vectors. The 
projections A and A' of the polarization vectors is 
conserved for S'-wave scattering. The expression given 
holds for A = A' = 0. This gives the only non- vanishing 
S-wave projection from the exchange diagrams: 



T 1 

J PS-Ex 



4^ 



w^ + w 3 ^ 2 



(v ■ k) 



yM\M 3 



(38) 



As explained in the main text, this contribution is for- 
mally of NLO, however, for the resonances discussed 
it gives a negligible contribution. 
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